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Reconstruction of lossless compressed fingerprint 
and x-ray images using Haar matrix technique 

Reem Talib Miri, Hussain Ali Hussain 


Abstract — The objective of this research is to decompress an 
image using steps not by one step to get an image with dimension 
and size different than the original without lossless any part of 
image the number of steps to reconstruction an image depend on 
the image dimensions. If the image dimensions are 512x512 then 
it needs 9 steps to decompress. And if it was 8x8 then it will need 
only 3 steps as explained in paper.. 

Index Terms —scaling function, Haar matrix, Haar wavelet, 
Lossless compression. 


i. Introduction 

From an historical point of view, wavelet analysis is a new 
method, though its mathematical underpinnings date back to 
the work of Joseph Fourier in the nineteenth century. Fourier 
laid the foundations with his theories of frequency analysis, 
which proved to be enormously important and influential. 
The attention of researchers gradually turned from 
frequency-based analysis to scale-based analysis when it 

started to become clear that an approach measuring average 
fluctuations at different scales might prove less sensitive to 
noise. The first recorded mention of what we now call a 
“wavelet” seems to be in 1909,in a thesis by Alfred Haar. The 
concept of wavelets in its present theoretical form was first 
proposed by Jean Morlet and the team at the Marseille 
Theoretical Physics Center working under Alex Grossman in 
France. The methods of wavelet analysis have been 
developed mainly by Y.Meyer and his colleagues, who have 
ensured the methods’ dissemination. The main algorithm 
dates back to the work of Stephane Mallat in 1988. Since 
then, research on wavelets has become international. Such 
research is particularly active in the United States, where it is 
spearheaded by the work of scientists such as Ingrid 
Daubechies, Ronald Coifman, and Victor Wickerhauser. 
Barbara Burke Hubbard describes the birth, the history, and 
the seminal concepts in a very clear text. See “The World 
According to Wavelets/’ A.K.Peters, Wellesley, 1996.The 
wavelet domain is growing up very quickly. A lot of 
mathematical papers and practical trials are published every 
month [2-3-4]. 

II. Haar Wavelets Operational Matrix 

We shall begin with the more convenient way for 
representing Haar wavelets in computer and for xe[A,B] 
which was given by [1] and for this purpose we define the 
quantity M=2 J , where J is the maximal level of resolution and 
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divide the interval [A,B] into 2M subintervals of equal length; 
each subinterval has the length. 

Ac =(B - A)/2M. 


Two parameters are introduced the dilation parameter j for 
which j— 0,1,..., J and the translation parameter 
k=0,l, ...,m — 1 where m =2 J . The wavelets number i is 
identified as i=m+k+l the i th Haar wavelet is defined as: 

( 1, for x G 

Hx) = j -1., for* E [£ 2 (r'),£ 3 (r')L "" (1) 

V 0, elsewhere. 


Where: 

f L (0 = .4 +- 2kuAx f 
(0 = .4 + (2 k + t)pAx f 
(0 = .4 FSCfc-hl^Ax 


[\=M/m 

The case i= 1 corresponding to the scaling function 

Taking the collocation points as following 

Xi=A+(i-0.5)Ax, t=l,2,... ,2M — (3) 

And Haar vector H 2 m (a) are defined as: 

H 2 m (xd = [ hp,), h 2 (x,h 2M fe) ] T ™ (4) 

By letting A=0, B=1 (usually the Haar wavelets are defined 

for the interval xE 0,1] j and hence Ax = — in eq.(3]..We 

2 M 


define the 2M-square Haar matrix H 2 mx 2 m as: 

H 2M X2M = [H 2M A) H 2M A) ... H im fAp]. (5) 

Let J=2 where / is the maximal level of resolution and divide 
the interval [0,1] in to 8 subintervals of equal length; each 
subinterval has the length Ax = 1/ 3. 

Two parameters are introduced the dilation parameter j for 
which j=0,l,2 and the translation parameter k=0 when m=l 
and k=0,l when m-2 and k=0,l,2,3 when m-4. The wavelets 
number i is defined as i=m-\-k+l. 
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For i=2, m=l, k =0 


p=4, Ax = 1/3 

(2) =0, ^ 2 (2) =1/2, & (2) =1 



x e [° ' J ' 
xE t -1 ]' 

else where* 


For i=3, m=2, k =0 
ju=2. Ax = 1/3 
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(3) =0, & (3) =1/4, (3) =1/2 

1 x e 0,J], 

h 3 w = < i _! 

0 else where. 

For i=4, m=2, k=l 
p= 2 , Ax = 1/8 

5i (4) =1/2, ^2 (4) =3/4, & (4) =1 

1 Ell 

h 4 W = i _! xeg.l], 

0 Elsewhere. 

For i= 5, m=4, k=0 
(j= 2 , Ax = 1 /S 

(5) =0, ^ (5) =1/8, ^ (5) =1/4 
1 *6 0,1], 

h5W = l “I -T S [j • j] ' 

0 elsewhere. 

For i= 6, m=4, k= 1 
p=2, Ax = 1/3 

(6) =1/4, ^2 (6) =3/8,_^ 3 (6) =1/2 

1 

^'=1-1 «e-i- 

For i= 7, m=4, k=2 
(j=l, Ax = 1/S 

(7) =1/2, (7) =5/8, (7) =6/8=3/4 

1 xe t;]' 

= i -1 xeg,g ( 

0 elsewhere. 

For i= 8, m= 4, k=3 
(j=l, Ax = 1/S 

(8) =3/4, (8) =7/8, (8) =1 

i » E - a - 

= i -1 xeg.l], 

0 efse wJwrre. 

Taking the collocation points as following: 
xj=l/16 , x 2 =3/16, x 3 =5/16 , x 4 -7/16, 


x 5 =9/16, x 6 =ll/16, x 7 =13/16, x 8 =15/6 
Haar vector H 8 (x,) are defined as: 

H 8 (x;) = [ h I (x,), h 2 (x,), h 3 (x,), h 4 (x,), h 5 (x,), h 6 (x,-), h 7 (x,), 
h 8 (x,) ] T 

For the first Haar vector 

H 8 (-t) = [ hi (A, h 2 (~), h 3 (~), h 4 A, h 5 A, h 6 A, h 7 A 

lb lb lb lb lb lb lb lb 

, h 8 A ] T 

1 b 

Hg A = [1,1,1,0,1,0,0,0] T 

1 b 

For the second Haar vector 

H 8 A = [ h,A, h 2 A, h 3 A, h 4 A, h 5 A, h 6 A, h 7 A 

Id Id Id Id Id Id Id Id 

, h*(~) ] T 

1 b 

h 8 A = [ i, i, i, o, -i, o, o, o ] T 

1 b 


For the third Haar vector 


H 8 (~) = [ h|(~), h 2 (~), h 3 A, h 4 ( 77 ), h 5 A, h 6 ( 77 ), h ? A 

. ] T 


16 


16 


16 


16 


16 


H 8 (~) = [ 1 , 1 , -1,0,0,1,0,0 ] T 

For the fourth Haar vector 

H 8 A) = [ hA, hA), hA), hA), h 5 A), hA), hA) 

lb lb lb lb lb lb lb lb 

, h 8 e> ] T 

16 

hA) = [ 1 , 1 ,-l ,0,0,-l ,0,0] T 

1 b 

For the fifth Haar vector 

H 8 A = [ hA), hA), hA), h 4 A, h 5 A, hA), hA) 

lb lb lb lb lb lb lb lb 

, hA) ] T 

1 b 

hA) = [ 1 , -1,0,1,0,0,1,0 ] T 

16 

For the sixth Haar vector 

h 8 A=[ h A), h 2 A, h 3 A, hA, h 5 A, h 6 A, h 7 A 

lb lb lb lb lb lb lb lb 

, h 8 A ] T 

1 b 

h 8 A = [ 1 ,- 1 , 0,1 ,o,o,-i ,of 

1 b 

For the seventh Haar vector 


hA)=[ hi A, h 2 A, h 3 A, h 4 A, h 5 A, h 6 A, h 7 A 

lb lb lb lb lb lb lb lb 

, h 8 A ] T 

1 b 

h 8 A = [ 1 ,- 1 , 0 ,- 1 , 0 , 0 , 0,1 ] T 

1 b 


For the eighth Haar vector 


h 8 §=[ h 4 >, h 2 §, h 4 >, h 4 g), h 5 §, h 4 >, h 7 §) 

15 T 

• *«<»> ] 


H 8 § = [ 1 ,- 1 , 0 ,- 1 , 0 , 0 , 0,-1 ] T 


We define the 8 -square Haar matrix H 8x8 as: 

h 8 *b = [ H 8 h 8 A h 8 A H 8 (— 3 h 8 A h 8 A 

H 8 § h 4 >] 



1 1 

1 1 

1 1 

0 0 

1 -1 
0 0 

0 0 

0 0 


1111 1 I 1 

1 1 - 1-1 -1 -1 
- 1-10000 
0 0 1 1 - 1-1 
0 0 0 0 0 0 
1-1 00 0 0 

00 1-1 0 0 

0 0 0 0 1 - 1-1 


The Haar matrix is the multiplication of the following three 
matrices. H 8x8 = h 3 . h 2 hi 

Where hi, h 2 , h 2 are: 
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III. STEPS FOR RECONSTRUCTION GRAYSCALE IMAGE 


Figure (2) represents the matrix C 
First step for reconstruction 


We will explain process on 8 x 8 image matrix. The process 
can be generalized to 2 J x2 J matrix, where is j natural number. 
If there are 2 J elements in row and column in matrix, then the 
process of row and column will consist of j steps. In our case j 
=3. Let A is an image matrix given below by using the 
command (imread image) in MATLAB. 



Figure (1) represents the matrix A 
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And C is an image compression matrix 


X 

The product of (h 3 ) xCxh 3 = RSi is equivalent to the 
following procedures. For every columns of C we calculate 
the addition of, and subtraction between, the first two 
elements. And replace the original numbers with the above 
results. We apply the same procedure to all other columns to 
get matrix Bi.Then for every rows of Bx we calculate the 
addition of, and subtraction between, the first two elements. 
And replace the original numbers with the above results. We 
apply the same procedure to all other rows to get matrix RSi. 
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Figure (3) represents the matrix RSi 
Second step for reconstruction 
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The product of ( Pt 2 ) xRS{Xh 2 = RS 2 is equivalent to the 
following procedures. For every columns of RSi we calculate 
the addition of the first and third elements and their 
subtraction. We also calculate the addition of the second and 
fourth elements and their subtraction. And replace the original 
numbers with the above results. We apply the same procedure 
to all other columns to get matrix B 2 .Then for every rows of 
B 2 we calculate the addition of the first and third elements and 
their subtraction. We also calculate the addition of the second 
and fourth elements and their subtraction. We repeat the same 
procedure for all rows to obtain matrix RS 2 . 
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Figure (4) represents the matrix RS 2 


Third step for reconstruction 

The product of (/q) x RS 2 x h\ = RS 3 = A is equivalent to the 
following procedures. For every columns of RS 2 we calculate 
the addition and subtraction of, first and fifth elements, 
Second and sixth elements, third and seventh elements, fourth 
and eighth elements. To get new matrix B 3 Then for every 
rows of B 3 we calculate the addition and subtraction of first 
and fifth elements, second and sixth elements, third and 
seventh elements, fourth and eighth elements. To get matrix 
RS 3 = A. 


RS 3 = A= 
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Figure (5) represents the matrix RS 3 equal to A 

We considered a small matrix because it is difficult to show 
the results on the large size of matrices like a matrix of size 
512x512. Now we use fingerprint or x-ray image of size 
512x512 =2 9 x 2 9 . First we will extend the inverse Haar 
matrix up to the size of the fingerprint or x-ray image. 
Later, we will apply the inverse Haar matrix to the image. The 
whole procedure will be carried out with the help of 
MATLAB. In our case the image matrix is of the dimension 
512x512 =2 9 x 2 9 , so inverse Haar matrix method 
compressed the image of fingerprint or x-ray in ninth steps 

o o 

the first step compresses to 256x256 = 2 x 2 the second 
step compresses to 128x128 = 2 7 x 2 7 the third step 
compresses to 64x64 = 2 6 x 2 6 the fourth step compresses to 
32x32 = 2 5 x 2 5 the fifth step compresses to 16x16 = 2 4 x 2 4 
the sixth step compresses to 8x8 = 2 x 2 the seventh step 
compresses to 4x4 = 2x2 the eighth step compresses to 2x2 
= 2 1 x 2 1 the ninth step compresses to lxl = 2° x 2° then we 
can sent the image by mail in internet or load it in website to 
decompressed the image after downloaded from the 
internet. We will used Haar matrix to reconstruction image 
in ninth steps the first step reconstructs to 2x2 = 2 1 x 2 1 the 
second step reconstructs to 4x4 = 2 x 2 the third step 
reconstructs to 8x8 = 2 x 2 the fourth step reconstructs to 
16x16 = 2 4 x 2 4 the fifth step reconstructs to 32x32 = 2 5 x 
2 5 the sixth step reconstructs to 64x64 = 2 6 x 2 6 the seventh 
step reconstructs to 128x128 = 2 7 x 2 7 the eighth step 

o o 

reconstructs to 256x256 = 2 x 2 the ninth step reconstructs 
to 512x512 =2 9 x 2 9 and we have the image. 
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